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The radiative return method - a short theory review 

Henryk Czyz''* 

'^Institute of Physics, University of Silesia, PL-40007 Katowice, Poland. 

A short review of the status of the theoretical developments concerning the radiative return method is presented. 
The emphasis is on the construction of the PHOKHARA Monte Carlo event generator and its tests. It is advocated 
that the radiative return method provides not only with the hadronic cross section extraction competitive with 
and complementary to the scan method, but also that it is a powerful tool in detailed studies of the hadron 
interactions. 



1. INTRODUCTION 

The hadronic cross section measurement is cru- 
cial for the accurate evaluation of the hadronic 
contributions to the muon anomalous magnetic 
moment (a^) [J and running of the electromag- 
netic coupling uqed Q]' The traditional way of 
measuring of the hadronic cross section via the 
energy scan has one disadvantage - it needs dedi- 
cated experiments. An alternative way, the radia- 
tive return method, was proposed in even if 
the radiative process was investigated earlier 0]. 
This method, described in the next section, al- 
lows for a simultaneous extraction of the hadronic 
cross section from the nominal energy of the ex- 
periment down to the production threshold, and 
importantly can profit from the data of all high 
luminosity meson factories. 



2. THE HADRONIC CROSS SECTION 
VIA THE RADIATIVE RETURN 
METHOD 

The radiative return method relies on an ob- 
servation that the cross section of the reaction 
e^e~ — > hadrons+photons, with photons emitted 
from the initial leptons, factorizes into a function 
H, fully calculable within QED, and the cross 
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section of the reaction 



H ■ da{e+e- 



— !■ hadrons 

-7's)(s,Q2) = 
hadrons) (Q2) , (1) 



where is the invariant mass of the hadronic 
system. Thus from the measured differential, 
in , cross section of the reaction e^e~ — s- 
hadrons -I- photons one can evaluate cr(e^e^ — > 
hadrons) once the function H is known. As evi- 
dent from the Eq. , the radiative return method 
allows for the extraction of the hadronic cross 
section from the production energy threshold of 
a given hadronic channel almost to the nominal 
energy of a given experiment (-\/s). The smaller 
cross section of the radiative process as compared 
to the process without photons emission has to 
be compensated by higher luminosities. That re- 
quirement is met by meson factories (DAPHNE, 
BaBar, BELLE). All of them were built for other 
purposes then the hadronic cross section mea- 
surements, but their huge luminosities provide 
with data samples large enough for very accu- 
rate measurements of interesting hadronic chan- 
nels and/or give an information on rare chan- 
nels, which were not accessible in scan experi- 
ments. Two representative examples of such mea- 
surements are the very accurate pion form fac- 
tor extraction by KLOE collaboration |5j and 
a{e^e~' — > Stt) extraction by BaBar collaboration 
ini, where it was shown that the old DM2 scan 
data were wrong at high values of . An exten- 
sive review of the recent results of both collabora- 
tions concerning the radiative return is presented 
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in [7]. If one likes to use the formula Eq.(Q in 
a realistic experimental situation, where sophisti- 
cated event selections are used, one needs a Monte 
Carlo event generator of the measured process. 
To meet that requirement the PHOKHARA |S] 
event generator was constructed. 

3. THE PHOKHARA MONTE CARLO 
EVENT GENERATOR AND ITS 
TESTS 

The construction of the PHOKHARA event 
generator started from the EVA generator 
where structure function method was used to 
model multi-photon emission. The physical ac- 
curacy of the program was however far from the 
demanding experimental accuracy of the KLOE 
pion form factor measurement and in a series of 
papers that high expectations were met. The first 
distributed PHOKHARA version relied on the 
one loop initial state radiative corrections calcu- 
lated in 11 and the two hard photon emission 
was simulated using exact matrix element written 
within helicity amplitudes method. That version 
was designed to run with tagged photon configu- 
rations and the radiative corrections necessary for 
a photon emitted at small angles were calculated 
afterwords in |12| and implemented into the event 
generator in |13) . The important issue of the fi- 
nal state emission, which will be discussed in de- 
tails in the next section, was addressed in and 
subsequently in [T5) . while in aspects specific for 
0- factory DAPHNE in |TB]. In parallel the gen- 
erator was being extended to allow for the gen- 
eration of more hadronic channels and now it al- 
lows for generation of tt+tt^, K^K^ , K'^K'^, pp, 
nn, 7r+7r~7r°, 27r+27r~, 7r+7r^27r*' hadronic states 
and fi^ fi^ . The nucleons final states were dis- 
cussed in 22], while the three pion current was 
modeled and implemented into the PHOKHARA 
event generator in |18j . 

All that allowed for building of the state-of-the- 
art event generator. The proper implementation 
of the radiative corrections as well as the hadronic 
currents is guarantied by extensive tests of the 
generator discussed below. 

At each step of the generator development, the 
newly implemented matrix element calculated in 



the program using the helicity amplitude method, 
squared and summed over all available helicities is 
compared with the square of the matrix element 
summed over polarizations calculated by tradi- 
tional trace method. All numerical calculations in 
PHOKHARA are performed using double preci- 
sion, but in some cases, mainly for double photon 
emission it was necessary to use the quadrupole 
precision for the matrix element evaluation cal- 
culated analytically by the trace method. It was 
caused by numerical cancellations up to ten sig- 
nificant digits occurring between various terms. 

Another type of tests concern the process of the 
generation. The initial state emission of one pho- 
ton with one-loop radiative corrections and two 
hard real photon emission from initial states were 
compared (13) separately with existing analytical 
results for fully inclusive phase space configura- 
tions '19' . Both results were in perfect agreement 
up to the numerical precision of the tests limited 
by the Monte Carlo statistics. The relative dif- 
ference between the numerical and analytical re- 
sults was a few times 10^"*, that was well within 
the statistical Monte Carlo error bars. That ac- 
curacy is usually called technical precision of the 
Monte Carlo generator and I will use that name 
hereafter. That tests are repeated for each newly 
added hadronic channel, even if the program uses 
the same building blocks for every channel, to 
avoid possible bugs in the implementation. 

Similar tests were performed for the reaction 
e+e^ — !■ 7r+7r~7(7) with one photon emitted from 
the initial state and one emitted from the final 
pions IE]- In that case the analytical results of 
were used together with the analytical re- 
sults obtained in [TJ to test the Monte Carlo 
generation. Again the same technical precision 
was achieved and the tests were repeated for the 
charged kaons in the final states. For the reac- 
tion e+e~ — !■ fi'^ fj,~ the configurations with 
one photon emitted from the initial states and 
one from final muons were tested in with the 
same technical precision. 

It would be useful to make extensive compar- 
isons with independent Monte Carlo generators, 
however the only existing Monte Carlo code meet- 
ing the accuracy requirements is the KKMC [5T] , 
which is limited to muons in the final state as 
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far as its accurate matrix element is concerned. 
It means that one can tests the initial state emis- 
sion and in fact detailed tests were performed [221 
leading to an excellent agreement of the non ex- 
ponentiated matrix elements of the virtual correc- 
tions to single photon emission (a relative differ- 
ence of a few times 10~^ was found). The higher 
order effects, that can be seen as a difference be- 
tween the exponentiated and the non exponenti- 
ated matrix elements reach at most 2 per mile 
with the exception of the region of the invari- 
ant mass of the hadronic system very close to the 
nominal energy of the experiment. That region, 
where soft multi-photon emission play an impor- 
tant role and thus the exponentiation is neces- 
sary, is however of no interest to the radiative re- 
turn method. All that results agree very well with 
the estimated previously in by means of the 
structure function approach, PHOKHARA phys- 
ical precision of 0.5%, attributed to the lack of the 
higher order effects in the ISR matrix element. 

4. THE FINAL STATE EMISSION 

The final state emission (FSR) forms a poten- 
tial problem for the application of the radiative 
return method and it has to be studied care- 
fully to be sure that the required accuracy of 
the description is met. First of all one has to 
have in mind that the situation at B- factories is 
completely different from the one of the (f)- fac- 
tory DAPHNE. In the former case the region of 
hadronic masses, which is of physical interests, 
mainly below 4 GeV, lays far from the nominal 
energy of the experiments, thus an emission of a 
hard photon is required to reach it. As a result 
the typical kinematic configuration of an event 
consists of a photon emitted in one direction and 
hadrons going opposite to it. That provides a nat- 
ural suppression of the FSR contributions, which 
are large for photons emitted parallel to the di- 
rection of a charged hadron in the final state, and 
makes the measurement of the hadronic cross sec- 
tion easier. For the 0- factory, where the physi- 
cally interesting region is not far from the nominal 
energy of the experiment, that natural separation 
between the emitted photon and the hadrons docs 
not exists and one has to suppress FSR by an ap- 



propriate event selection. In that case one has 
to control the uncertainty due to the model de- 
pendence of the final state emission. That is a 
challenge, as the models were not tested with the 
adequate precision prior to the DAPHNE results. 
I will discuss that problem on the basis of the 
e^e~ — > 7r^7r~7(7) process, where the accuracy 
requirements are the most demanding. For that 
process the solution of the problem was first pro- 
posed in and further elaborated in ^3]. One 
can imagine a similar solution for other hadronic 
final states, but that kind of analysis was never 
performed. 

The main tool in the tests of the model(s) of 
the photon emission from the final pions is the 
charge asymmetry. For one real photon emission 
the two-pion state is produced in C=-l and with 
an odd orbital angular momentum for the real 
photon emitted from the initial state and in C=l 
and with an even orbital angular momentum for 
the real photon emitted from the final state. As 
a result, the initial- final state interference is odd 
under 7r+ 7r~ interchange and it integrates to 
zero for charge blind event selections. In the same 
time it is the only source of the charge asymmetry 
and as such allow for tests of the models of the 
final state emission. The charge asymmetry de- 
pends on the invariant mass of the two-pion sys- 
tem and that allow for deeper insight into details 
of the tested model(s). In short, the tests should 
be done in the following way: First one compares 
the experimental data for the asymmetry with the 
Monte Carlo where the tested model was imple- 
mented. That has to be performed for an event 
selection which enhance the FSR as compared to 
the ISR. Once the implemented model agrees with 
the data one chooses an event selection, which 
suppresses the FSR and performs the radiative 
cross section measurement. That guaranties that 
the ISR and the FSR contributions are separately 
well under control. For the case of untagged pho- 
tons a specific background, e~^e~ — > 7r"'"7r~e"'"e~, 
has to be also taken into account |23I24| as the 
final leptons are not vetoed. 

The reaction e"'"e^ tt^tt^j, with the pho- 
ton emitted from the pions, does contribute also 
to dispersion integrals for evaluation of and 
oiQED and in the former case its theoretically es- 
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timated value is of the size of the theoretical 
uncertainty and thus numerically important. As 
its theoretical estimations are not reliable it has 
to be measured. The sketched program was suc- 
cessfully undertaken by KLOE and resulted in a 
sound extraction of the cr(e+e^ tt+tt^) 5 to- 
gether with the mentioned FSR photon correc- 
tions. 

Another source of complications for using of 
the radiative return method at DAPHNE are the 
radiative (j) decays. That problem was considered 
for the first time in |25j and it is discussed in more 
details in the next section. 

5. THE RADIATIVE RETURN AS A 
TOOL IN HADRONIC PHYSICS 

The reaction e'^e' — > ((> — > 7r+7r^7 produces 
the same final state as the one measured for the 
a{e~^e~ — > 7r"'"7r~) extraction. That contribution 
is sizable for energy close to the (f) mass and thus 
important for DAPHNE. As shown in IB', the 
charge asymmetry has large analyzing power and 
can provide with information allowing for distin- 
guishing between different models of the radiative 
(j) —^ TTTTj decay, even if that is impossible in the 
analysis of the differential (in Q^) cross section. 
Again by an appropriate event selection one can 
suppress those contributions or enhance them as 
for other sources of the FSR emission discussed 
in the previous section. That example shows that 
the radiative return method can be used not only 
for the hadronic cross section measurement, but 
also for getting detailed information about the 
models of hadronic interactions. It was partly 
exploited in the KLOE analysis however as 
the charge asymmetries were not used in the fits, 
the collected data contain more information on 
the tested models then actually was used. 

An extensive analysis of the FSR contributions 
is extremely important especially in the ongo- 
ing KLOE analysis, both for tagged and un- 
tagged j2Hl photon(s) as till now the experimen- 
tal information on the pion-photon interactions is 
far from being satisfactory and it is not clear if 
the model used currently in PHOKHARA (sQED 
+ vector dominance + radiative (j) decays) de- 
scribes the FSR with the adequate precision in 



the threshold region, where other contributions 
might be important |29| . 

Another example of the power of the radiative 
return method in the hadronic models tests is the 
separation of the magnetic and the electric nu- 
cleon form factors. The method, which was pro- 
posed in ;l7i. was used by BaBar collaboration 
30 for separation of the proton form factors. The 
obtained results show clearly its competitiveness. 

6. THE SUMMARY AND NEAR FU- 
TURE DEVELOPMENTS 

A short description of the theoretical status 
of the radiative return method was presented 
showing its competitiveness in precise measure- 
ments of the hadronic cross section and stud- 
ies of the hadronic interactions. Many inter- 
esting problems, for example a proper model- 
ing of the hadronic current of multi-meson final 
states observed at BaBar 7 , the FSR simula- 
tion for more than the two-pions final states, the 
modeling of the narrow resonance contributions 
and many others not mentioned in this paper 
still await for detailed theoretical investigations. 
One of that problems, which will be addressed 
in the near future by the group working on the 
PHOKHARA event generator developments and 
updates, is the improvement of the theoretical 
description of the An hadronic current Ex- 
ploiting isospin symmetry and all available exper- 
imental data one comes to the predictions for the 
CT(e+e^ 27r*'7r+7r^) (central dashed line) shown 
in FigH The lower and upper dashed lines show 
the error bars of the model predictions. 

Acknowledgments: The author thanks the 
organizers of the The International Workshop 
"e-l-e- collisions from cj) to for hospitality and 
stimulating atmosphere of the whole event. 

REFERENCES 

1. M. Passera, these proceedings. 

2. F. Jegerlehner, these proceedings. 

3. Min-Shih Chen and P. M. Zerwas, Phys. Rev. 
D 11 (1975) 58. 

4. V. N. Baier and V. S. Fadin, Phys. Lett. 27B 
(1968) 223. 

5. A. Aloisio et ai, [KLOE CoUabora- 



5 



20 



10 











1 

' ^ \ 


1 

CMl 


JD - -- 
JD --- 
)-2 - --- 
G2 -- 
lei 


















vmi 


















W t 








— ' 


1 




1 







0.6 0.8 1 1.2 1.4 1.6 1.8 2 22 2.4 
£c,». IGeV] 



Figure 1. The predictions for the a{e^e~ 
27r°7r+7r~) (central dashed hne). The lower and 
upper dashed lines show the error bars of the 
model predictions. 



tion], Phys. Lett. B 606 (2005) 12, 
'hep-ex/0407048 . 

6. B. Aubert et al. [BaBar Collabora- 
tion] Phys. Rev. D 70 (2004) 072004 
'hep-ph/0408078 .. 

7. A. Denig, these proceedings 

8. G. Rodrigo, H. Czyz, J.H. Kiihn and 
M. Szopa, Eur. Phys. J. C 24 (2002) 71 
hep-ph/01 12184. 

9. S. Binncr, J. H. Kiihn and K. Melnikov, Phys. 
Lett. B 459 (1999) 279 hep-ph/9902399 . 

10. H. Czyz and J. H. Kiihn, Eur. Phys. J. C 18 
(2001) 497 hep-ph/0008262 . 

11. G. Rodrigo, A. Gehrmann-Dc Ridder, 
M. Guilleaume and J. H. Kiihn, Eur. Phys. 
J. C 22 (2001) 81 hep-ph/0106132 . 

12. J. H. Kiihn and G. Rodrigo, Eur. Phys. J. C 
25 (2002) 215 hep-ph/0204283 . 

13. H. Czyz, A. Grzelihska, J. H. Kiihn and 
G. Rodrigo, E ur. Phys. J. C 27 (2003) 563 
|hep-ph/ fe2225j . 

14. H. Czyz, A. Grzelihska, J. H. Kiihn and 
G. Rodrigo, Eur. Phys. J. C 33 (2004) 333 
[hep-ph/0308312 . 



15. H. Czyz, A. Grzelihska, J. H. Kiihn and 
G. Ro drigo, Eur. Phys. J. C 39 (2005) 411 
|hep^p h/0404078l. 

16. H. Czyz, A. Grzelihska and J. H. Kiihn Phys. 
Lett. B 611 (2005) 116 hep-ph/0412239 . 

17. H. Czyz, J. H. Kiihn, E. Nowak and G. Ro- 
drigo, Eur. Phys. J. C 35 (2004) 527 
hep-ph /04030 62 . 

18. H. Czy z, A. Gr zehhska, J. H. Kiih and G. Ro- 
drigo, hep-ph/0512180 

19. F. A. Berends, W. L. van Neerven and 
G. J. Burgers, Nucl. Phys. B 297 (1988) 429; 
Erratum, ibid. B 304 (1988) 921. 

20. J. S. Schwingcr, Particles, Sources, and 
Fields. Vol. 3, Redwood City, USA: Addison- 
Wesley (1989) p. 99. 

21. S. Jadach, B. F. L. Ward and Z. Was, 
Comput. Phys. Commun. 130 (2000) 260 
hep-ph/9912214 . 

22. S. Jadach, B. F. L. Ward and S. A. Yost, 
Phys. Rev. D 73 (2006) 073001 
lhep-ph/0602197 ; S. A. Yost, S. Jadach 
and B. F. L. Ward, Acta Phys. Polon. B 36 
(2005) 2379 hep-ph/0505065 ; S. Jadach, 
Acta Phys. Polon. B 36 (2005) 2387 
|hep-ph/0506180 ; C. Glosser, S. Jadach, 
B. F. L. Ward and S. A. Yost, Phys. Lett. B 
605 (2005) 123 hep-ph/0406298 . 

23. A. Hoefer, J. Gluza and F. Jegerlehner, 
Eur. Ph ys. J. C24 (2002) 51, 
"hep-ph/0107154 . 

24. H. Czyz and E. Nowak-Kubat, Phys. Lett. B 
634 (2006) 493 hep-ph/0601169 . 

25. K. Melnikov, F. Nguyen, B. Valeriani and 
G. Venanzoni, Phys. Lett. B 477 (2000) 114 
hep-ph/0001064 . 

26. F. Ambrosino et al. [KLOE Collabora- 
tion], Phys. Lett. B 634 (2006) 148 
jhep-ex/0511031 . 

27. D. Leone, these proceedings 

28. F. Nguyen, these proceedings 

29. G. Pancheri, O. Shekhovtsova and G. Venan- 
zoni, hep-ph/0605244, 

30. B. Aubert et al. [BABAR Collabora- 
tion], Phys. Rev. D 73 (2006) 012005 
arXiv:hep-ex/0512023 . 

31. H. Czyz, J. H. Kiihn and A. Wapienik, in 
preparation. 



